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The interference model has been widely used and studied in block 
experiments where the treatment for a particular plot has effects on 
its neighbor plots. In this paper, we study optimal circular designs 
for the proportional interference model, in which the neighbor effects 
of a treatment are proportional to its direct effect. Kiefer’s equiv¬ 
alence theorems for estimating both the direct and total treatment 
effects are developed with respect to the criteria of A, D, E and T. 

Parallel studies are carried out for the undirectional model, where 
the neighbor effects do not depend on whether they are from the 
left or right. Moreover, the connection between optimal designs for 
the directional and undiretional models is built. Importantly, one can 
easily develop a computer program for finding optimal designs based 
on these theorems. 

1. Introduction. In many agricultural experiments, the treatment as¬ 
signed to a particular plot could also have effects on its neighbor plots. This 
is well recognized in literature. See Draper and Guttman (1980), Kemp- 
ton (1982), Besag and Kempton (1986), Langton (1990), Gill (1993) and 
Goldringer, Brabant and Kempton (1994), for examples. To adjust the biases 
caused by these neighbor effects, the interference model is widely adopted. 
In a block design with n blocks of size k and t treatments, the response, 
yuij, observed from the jth plot of block i is decomposed into the following 
items: 

(1) Udij = UP Pi + Td(ij) + 7ci(i J_l) + Pd{i,j+l) + ) 

where the subscript d{i,j) denotes the treatment assigned to the jth. plot of 
block i by the design d: {1,2,..., n} x {1,2,..., A:} —>■ {1,2,..., t}. Further¬ 
more, /i is the general mean, (3i is the effect of block i, Td(ij) is the direct effect 
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of treatment d{i,j), Xd{i,j-i) is the neighbor effect of treatment d{i,j — 1) 
from the left and Pd{i,j+i) is the neighbor effect of treatment d{i,j + 1) from 
the right. At last, is the error term. 

Kunert and Martin (2000) studied optimal designs under model (1) for 
estimating the direct treatment effect when A: = 3 or 4. The latter was ex¬ 
tended to 5 < A; < f by Kunert and Mersmann (2011). Zheng (2015) recently 
provided a unified framework in deriving optimal designs for general val¬ 
ues of k and t, with an arbitrary structure of the within-block covariance 
matrix. On the other hand, Bailey and Druilhet (2004) studied the optimal 
designs under the same model, however, for estimating the total treatment 
effect which is the summation of the direct and neighbor effects. This line 
of research was extended by Ai, Ge and Chan (2007), Ai, Yu and He (2009) 
and Druilhet and Tinsson (2012). See also Gill (1993), Druilhet (1999), Fil- 
ipiak and Markiewicz (2003, 2005, 2007) and Filipiak (2012) among others 
for relevant works on optimal designs. 

In this paper, we shall consider the proportional interference model, where 
the neighbor effects are proportional to the direct treatment effect, that is, 
li = AiTi and pi = A 2 Tj, 1 < A < A, for unknown constants Ai and A 2 . This is 
reasonable for many applications since an effective treatment typically has 
large impacts on its neighbor plots. In fact. Draper and Guttman (1980) 
has proposed such model with Ai = A 2 . A model with this restriction is 
said to be undirectional; otherwise, it is directional. Yet, there is no lit¬ 
erature on optimal designs under either of these two models according to 
the best knowledge of the authors. Meanwhile, optimal crossover designs 
under a similar proportional model have been studied by Kempton, Ferris 
and David (2001), Bailey and Kunert (2006), Bose and Stufken (2007) and 
Zheng (2013a). By their enlightenment, the nonlinear terms XiTi and A 2 Tj 
in the proportional interference model can be handled in the same fashion. 
We are interested in finding the optimal designs for estimating the direct 
and total treatment effects, respectively, under either of the directional and 
undirectional models. 

Let Yd be the vector of responses organized block by block. Now we can 
write the proportional interference model as follows: 

(2) Yd = ^nkP + U(3 + {Td -|- XiLd + X2Rd)T + £; 

where Ink represents a vector of ones with length nk, [3 = (/3i,...,/!„)', 
r = (ri,..., Tt)' and U = In ® '^k- Here, In is the identity matrix of order 
n; (D denotes the Kronecker product and ' means the transposition. Also, 
Td, Ld and Rd represent the design matrices for the direct, left neighbor 
and right neighbor effects, respectively. Throughout the paper, we con¬ 
sider circular designs, for which d{i,0) = d{i,k) and d{i,k -|- 1) = d{i,\), 
1 <i <n. Hence we have Ld = {In <S> II)Td and Rd = {In ® II')Td, where 
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H = (IIi=j+i(mod k))i<ij<k with the indicator function I. For the random er¬ 
ror term e, we assume that E(e) = 0 and Var(e) = /„ (8) S, where S is an 
arbitrary k x k positive definite within-block covariance matrix. 

The rest of the paper is organized as follows. Sections 2 and 3 investigate 
the optimal designs for estimating the direct and total treatment effects, 
respectively, under the proportional interference model. Kiefer’s equivalence 
theorems are given with respect to A, D, E and T criteria therein. Sec¬ 
tion 4 carries out parallel studies for the undirectional model. Moreover, the 
connection between optimal designs for the two models is built. Section 5 
illustrates these theorems through several examples. Section 6 concludes the 
paper with some discussions. 


2. Optimal designs for direct treatment effect. For any matrix G, de¬ 
fine G~ as a generalized inverse of G and the projection operator pr-*- G = 
I - G{G'G)-G'. Let K = (4 ® , 4 = (4 ® = [In ® 

and Rd = (4 ® The Fisher information matrix for the 

direct treatment effect r under model (2) is 

Cd{T) = (4 + AiLrf -|- X2Rdy P'^"^(U\LdT\RdT){Td -|- XiLd + \2Rd)- 

For notational convenience, let = xTd + yLd + zRd for any values of 

x,y and z. By setting Aq = 1, we have 

Cdir) = ^i^{U\LdT\RdT)Mi^XuX2 

= Pr-^(^)Tfl,Ai,A2 - Tf{,Ai,A2 pr-^{U){LdT\RdT) 

( 3 ) x[{LdT\RdTy pv-^{U){LdT\RdT)]~ {LdT\RdTy pv-^{U)Mi^Xi,\2 


Ad 


2 2 

EE K^jCdij A^{t 

j=0 


XiCdilT 




2 = 0 


where Gdij = G[{In ® B)Gj, 0 < i,j < 2, with Go = 4, Gi = Ld, G 2 = Rd 
and B = In particular, if S is a matrix of 

type-H, that is, S = 4 + V^'k + ^kf]' with a vector y of length k, we have 
B = pr-’-(lfc) := Bk [Kushner (1997)]. Examples of type-H matrices include 
the identity matrices and completely symmetric matrices. 

One major objective of design theorists is to hnd a design with maximum 
information matrix. Following Kiefer (1975), we shall try to find the designs 
which maximize 4>(Grf(T)), where $ satishes the following three conditions: 

(C.l) is concave; 
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(C.2) ^{S'CS) = for any permutation matrix 5; 

(C.3) is nondecreasing in the scalar b> 0. 

Note that Cdij) depends on the true value of r itself, and thus the choice 
of optimal designs. Following Kempton, Ferris and David (2001), Bailey and 
Kunert (2006) and Zheng (2013a), we adopt the Bayesian type criterion 

(4) Ud) = I ^Cd{T))g{T) d(T) = EgmCdir))), 

where g is the prior distribution of r and is assumed to be exchangeable 
throughout the paper. A design is said to be optimal if and only if it achieves 
the maximum of 4>g{d) among all designs for given gf, <h, Ai and A 2 . Further¬ 
more, if a design maximizes 4>g{d) for any <h, it is also said to be universally 
optimal. 

In this paper we consider four popular criteria for finding optimal designs. 
For a t X t matrix C with eigenvalues 0 = oq < oi < 02 < ■ • ■ < o-t-i-, define 
the criterion functions as 

^A{C) = {t-l)C^a-^ 

Vi=i 
n-i \ V(i-i) 

‘i>d(c)= 

^e{C) = oi, 

t-i 

$r(C') = (t-l)-'j;a,. 

i=l 

A design is said to be Ag-optimal if it maximizes 4>g{d) with in (4). 

The Vg-, Eg- and 7^-optimality of a design are similarly defined. 

Let Eln^k,t denote the set of all possible block designs with n blocks of size 
k and t treatments. A design in could be considered as a result of 

selecting n elements from the set, 5, of all possible block sequences with 
replacement. For each s a S, we define the sequence proportion ps = ngjn, 
where Ug is the number of replications of s in the design. For given n, a design 
is determined by the measure ^ = {ps,s G S). ps > 0, then s is a supporting 
sequence of In approximate design theory, we search for optimal measures 
in the space of T* = {C-YlsesP^ ~ such a measure happens 

to fall within the subset Vn = {^ G V :n^ is a vector of integers}, then we 
derive an exact design which is optimal among Eln,k,t- 

Let Csij be the matrix Cdij when the design d is degenerated to a single 
sequence s for 0 < i,j < 2. Then we have Cdij = ^By equation 
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(3), we have 


Cd{T) =nC^{T), 


2 2 


(5) 


XiXjC^ij A^{t 

i=0 j=0 


'Sljl 


^iij — 'y ^j PsCs 
s^S 


0=0 


i=0 


Note that C^(t) is independent of n. Here we call C^{t) the information 
matrix of the measure Furthermore, by noting that the four criterion 
functions satisfy ^{nC) = n^{C), we have 


( 6 ) 


(t)g{d)=n(t)g{i), 

(t^giO = J ^{CdT))9{r)d{T). 


Equation (6) indicates that the number of blocks n is irrelevant to the search 
of approximate optimal designs. In the sequel, we shall focus on Ending the 
optimal measures which maximize (pg{^) among V. 

Let O denote the set of all t\ permutation operators on {1,2,..., t}. For 
any a & 0 and s = {ti,... ,tk) with 1 < L < t, define as = {a{ti ),..., a{tk))- 
A measure is said to be symmetric if it is invariant under treatment rela¬ 
beling, that is, a^ = ^ for all u G O, where a^ = (Po-~is, s G 5). By adopting 
similar arguments to Corollary 1 in Zheng (2013a), we get the following 
result. 


Proposition 1. In approximate design theory, given any values of Xi 
and X 2 , and the exchangeable prior distribution g of t, for any measure f 
there exists a symmetric measure, say If, such that 

4>g{0 < 4>g{C)- 

Proposition 1 indicates that an optimal measure in the subclass of sym¬ 
metric measures is automatically optimal among V. The merit of such a 
result is that the form of the information matrix for a symmetric measure is 
usually feasible to be calculated explicitly. In fact there is a larger subclass of 
measures with the same convenience. We say a measure is pseudo symmetric 
if C^ij, 0 <i,j <2 are all completely symmetric. A symmetric measure is 
also pseudo symmetric [Kushner (1997)]. It is easy to verify that the column 
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and row sums of C^ij’s are all zero. Hence, for any pseudo symmetric mea¬ 
sure we have C^ij = c^ijBt/{t — 1), 0 < i,j < 2, where c^ij = tr(C'^jj). Now 
let £= (l,Ai,A 2 y, = (c^ij)o<i,j< 2 , 

( 7 ) g| = c^oo - (c^oi c^02) Qy ^ • 


Proposition 2. For a pseudo symmetric measure the information 
matrix C^{t) has eigenvalues o/O, {t — and {t — with mul¬ 

tiplicities of 1, 1 and t — 2, respectively. Moreover we have 


Proof. Due to 1[t = 0, we have Btr = r and r'C^ijr = c^i^Fr j (t — 1). 
In view of (5), we obtain 


2 2 


{t-C)C,(T) = Y£L XiXjC^ijBt - air'r) Vr', 

i=0 j=0 

/ 2 2 \ 

^ y ^ y XiC^i2 I 2 

\ i=o / 


a = 


^ 2=0 2=0 


Let {xi,... ,xt- 2 } be the orthogonal basis that is orthogonal to both It 
and T. Then {If, r, xi,..., Xi_ 2 } forms the eigenvectors of C^{t). The cor¬ 
responding eigenvalues are 0, {t — iYl‘i=oYl‘j=o ~ “ 

1)“^ AjAjC^fj with multiplicities of 1, 1 and t — 2, respectively. 

The proof is concluded in view of 


2 2 


i V^l — XiXjC^ij, 

i=0 j=0 


and (7). □ 


/ c^oi 

c^ii 

Vc«2i 


c?02\ 

C 512 

C 522 / 



C 511 

C 512 


V > 0 

C 522 / 


By Proposition 2, it is seen that 4>g{f,) = for any pseudo symmet¬ 

ric measure under the four criterion functions. Hence g is irrelevant to the 
determination of optimal pseudo symmetric measures for the four criteria. 


Lemma 1 . Except for measures with each supporting sequence consisting 
of only one treatment, we have c^a > 0 for i = 0,1,2. If det((5^) > 0, then 
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= det(V^)/det(Q^), where det(-) means the determinant of a matrix. Oth¬ 
erwise, = c^oo — c|oi/c 5 ii = c^oo — c|o 2 /c^ 22 - 

Proof. Note that B is nonnegative definite. So c^oo = 

> 0, where Tg is the matrix when d is degenerated to a single sequence 
s. If c^oo = 0) we have T^BTs = 0 and thus BTg = 0 for any supporting 
sequence s. It is known that Bx = 0 if and only if x is a multiple of 1^ 
[Kushner (1997)]. This is only possible when each supporting sequence re¬ 
peats the same treatment throughout the k plots. For c^n and 0 ^ 22 , we have 
the similar arguments. The rest of the lemma follows by straightforward 
calculations. □ 


From the proof of Lemma 1, we know Vg = 0 if each supporting sequence 
of f consists of only one treatment. There is no information gathered from 
such measures regarding r, and hence it is impossible to be optimal. In the 
subsequent arguments, we neglect such measures by default. 

Define the quadratic functions q^{x) = cgoo + 2cgoi^ + cgnx^ and xg = 
-cgoi/c^ii- Then we have q^{x^) = cgoo - c|oi/c 5 ii- Let Cgij = tr(C'sij), K = 
(Csij)o<i,j<2) Qs — (c and qs{x) = c^oo + 2csoia; + c^nx^. Clearly, 

Qi = Y.s&sPsQs and q^{x) = YlsesPsQsix) 


Theorem 1. In estimating r under model (2), a pseudo symmetric mea¬ 
sure f is optimal in the following cases. In each case, the Ps in ^ is positive 
only if s reaches the maximum therein. 


(i) //det(Qg) = 0, then f is Ag-optimal if and only if 

q^{x^)-^qs{x^) + it- 

//det(V^) > 0, then ^ is Ag-optimal if and only if 

r,det(Qg)/det(yg) + (t - 
ses detlQ^)/det(V^)-h (t — 2)(£'V^£)~^ 

(ii) If det{Q^) = 0, then f is Vg-optimal if and only if 


( 8 ) 


/ 1 g,(xg) t-2£'Vg£ \ 
\t — lq^{x^) t — lBV^l) 


//det(V^) > 0, then ^ is Vg-optimal if and only if 


(9) 


max 

sGS 


rg t-2eVg£\_ 
t-1^ t-1 £'V^£) 
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(iii) If det{Q^) = 0, then f is £g-optimal if and only if 


Qs{x^) 
max ——^ 
ses q^{x^) 


= 1 . 


//det(V^) > 0, then ^ is £g-optimal if and only if 


maxr, = 1. 

ses 

(iv) If det{Q^) = 0, then ^ is Tg-optimal if and only if 

qs{x^) + {t-2yVsI 
q^{x^) + {t- 2)£'V^e 

//det(V^) > 0, then ^ is Tg-optimal if and only if 

^ rsdet(V"^)/det((5g) + (i - 2)i'Vsi _ ^ 
det(Vg)/det(g5) + (i - 2)£'V^£ ~ 

Here r* = tr(T4y^“^) - tr(QsQ^^). 


Proof. Here we give only the proof for (ii) and the other three cases 
follow similarly. First we would like to show that 

(10) det(P^)/det(g^) < c^oo - c|oi/c^ii 

whenever det(gg) > 0. To see this, consider the following inequality: 

(11) J_/'c^02\( 

V'^ao c^iiy c^22V'^«i2y ^ ^ c^iiy 

The left (resp., right) hand side of (10) is the Schur complement of the left 
(resp., right) hand side of (11), and hence (10) follows by the nondecreasing 
property of Schur complement. 

By the definition of P^-optimality, Propositions 1 and 2, Lemma 1 and 
inequality (10), a pseudo symmetric measure I, with det(g^) = 0 is Vg- 
optimal if and only if 


( 12 ) 


hm V^[(l-'^)^ + ^^o]-^(0 

( 5->-0 h 


for any measure .^Oj where = log(( 7 ^(xg)) + {t — 2) log(^'Vg^). Here we 
used the fact that V^* = and hence ^'(Co) ~ '4’if.o)^ where is a sym¬ 
metric measure defined by Direct calculations show that 

(12) is equivalent to 

1 ggo(^g) t-21%^1 ^ 
t — 1 q^ix^) t — \ ~ 


(13) 
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By reducing to a degenerate measure which puts all mass on a single 
sequence s, we have 


(14) 


max 

sGS 


1 g,(xg) t-2£'VJ \ ^ 
t — 1 g^(x^) t — 1 £'V^£J ~ 


By letting = we have equality in (13) and hence 


(15) 


max 

sGS 


1 qs(x^) ^ 

t — 1 g^(x^) t — 1 £'V^£) ~ 


in view of g^{x) = XsG 5 Ps^s(a^)- Combining (14) and (15), we obtain (8). 

For a pseudo symmetric measure ^ with det(V^) > 0 and any measure ^Oj 
by the continuity of det(( 5 (i_ 5 )^_|_, 5 ^g) in 6, there exists a constant e > 0 such 
that det((5(i_5)^+5^o) ^ ^ ^ Hence ^ is D^-optimal if and 

only if 


(16) 


lim + 

5-i^O 6 


where (p{^) = log(det(Vg)/det(( 5 ^)) + {t — 2 ) log(^'Vg^). It is well known that 

/ N log(det(V’n_5)c 15f„)) — log(det(V£)) ^ , 

17 lim^^- ^ (1 -Bv-^ y y-i _ 3_ 

^ ' 5^-0 6 V ?o 4 / 


The same result holds for except that the number 3 in (17) is replaced 
with 2. By applying (17) to (16) we have 


H(HgoVg ^)-tr(Qgog 


(18) 


-i> 


f- 1 

Hence, for single sequences we have 

r-l 


f - 1 £'V^£ - ■ 


max 
s^S V 


/tr(F^F^ ) - tiiQsQ^ ) -2 £'Vs£ 


t-l 


t - 1 £'V^£ 


< 1 . 


By taking we have equality in (18). Also observe that conditioning 

on fixed the left-hand side of (18) is a linear function of the proportions 
in ^ 0 - Hence we have 


max 

sG.S 


tr(14F^ ^) - tr(Q^Q^ ^) t-2 £'Vs£\ 

T^i ^T^Jv[£)- ■ 


Then equation (9) follows. □ 


Remark 1. Theorem 1 neglected the pseudo symmetric measures with 
det(V^) = 0 and det((5^) > 0, and Theorem l(i)-(iii) neglected those with 
det(Qg) = 0 and g^{x^) = 0. However, all these measures yield g"^ = 0, and 
thns they cannot be optimal under Ag, Vg and £g criteria. Actually, r is 
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not estimable for any measure with = 0 and hence such measures should 
not be adopted [Pukelsheim (1993), Chapter 3]. Note also that det(Q^) = 0 
implies det(V^) = 0. Theorem 1 gives a comprehensive list of conditions to 
judge the optimality of a pseudo symmetric measure for estimating r. 

Remark 2. Since we also have = c^oo “ c|q 2 /c ^22 by Lemma 1, if the 
function q^{x) is replaced with c^qo + + 0 ^ 22 ^'^^ equivalent conditions 

for optimal pseudo symmetric measures with respect to the four criteria 
could be derived similarly. 

Remark 3. For the nonproportional model (1), the information matrix 
of a pseudo symmetric measure has t — 1 eigenvalues of q’^ and one of 0. The 
measure in Theorem l(iii) should also be universally optimal under model 
( 1 ). 

Remark 4. When there is only one neighbor effect, say left, we have 
det((5g) = 0 for all ^ £V. Theorem 1 reduces to equivalent conditions for 
the optimal crossover measures where the pre-period treatment is equal to 
the treatment in the last period for each subject. 

3. Optimal designs for total treatment effect. In this section we study 
optimal measures for estimating the total treatment effect, dehned hy 6 = 
(1 -|- Ai -|- A 2 )t. Bailey and Druilhet (2004) commented that the total treat¬ 
ment effect is more important when the experiment is aimed at finding a 
single treatment which is recommended for use in the whole held. 

When 1 -|- Ai -I- A 2 = 0, 0 takes the value of constant 0 regardless the value 
of r, and there is no need to carryout the experiment. In the following we 
assume 1 -|- Ai -I- A 2 7 ^ 0. By plugging t = 9/{1 + \i + A 2 ) into model (2), we 
have 

Fd = Infc/^ -|--|- (1 -|- Ai -|- A 2 ) -|- AiLrf + ^2Rd)9 + £■ 

The information matrix for 6 is 

Cd{9) = (1 -|- Ai -|- A 2 ) ^ 

^ ^ 1 ,Ai,A 2 PI"^(^l^-l,l+A 2 ,-A 2 ^l^-l,-Ai,l+Ai 6 *)Mi^Ai,A 2 - 

Here we used the equation pr-*- EF = pr-*- E for any nonsingular matrix F. 
Actually, it is seen that 1 -|- Ai -|- A 2 = 0 will yield inhnite Cd{9) for any d, 
which implies that the covariance matrix for 0 is zero. Our previous comment 
on this special case is justihed here. In the same way that we dehned C^{t) 
in Section 2, the information matrix of a measure ^ for 9 is given by C^{9) = 
n“^C'rf( 0 ), which is independent of n and can be expressed in a similar 
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fashion to equation (5). In the spirit of Proposition 1, we shall restrict our 
considerations to pseudo symmetric measures. 

To precede, we define £o = (“Ij 1 + -^ 2 , —A 2 )^ = (—1, —Ai, 1 + Ai)', Lq = 
(4,4) and Li = (£,4,4). Let = L[V^Li, and = 

£%£-£%LoQl-^L'oV^e. 

Proposition 3. For a pseudo symmetric measure 4 the information 
matrix C^{6) has eigenvalues of 0, (1 + Ai + A 2 )“^(i — and (1 + 

Ai + A 2 )~^(t — 1)~^£'V^£ with multiplicities of 1, 1 and t — 2, respectively. 
Moreover we have i< £'V^£. 

Proof. Denote Ad = (M_i,i+A 2 ,-A 2 ^'|Af-i-A i.i+Ai^*). Using 1(0 = 0 and 
C^ij = c^ijBt/{t - 1), we have 

C^{9) = n ^(1 + Ai + A 2 ) 

X {^LAi,A 2 Pr^(^)ALl,Ai,A2 

-Km,\2 Pr^(^)^rf[^dPr^(^)^rf]"^dPr^(^)ALi,Ai,A2} 

= (1 + Ai + X2)-\t - l)-^[£%£Bt - a{6'6)-^ee'], 

where a = ^LqV^£. Let {xi,..., xt- 2 } be the orthogonal basis which 

is orthogonal to both It and 9. Then {li,0,xi,... ,xt_ 2 } forms the eigen¬ 
vectors of C^{9). The corresponding eigenvalues are 0, (1 -|- Ai -|- A 2 )“^(t — 
and (1 -|- Ai -|- \ 2 )~‘^{t — 1)“^/Vg£ with multiplicities of 1, 1 and 
t — 2, respectively. The proof is completed in view of a > 0. □ 

Since Vg = 0 implies Vgq = 0, we neglect the measures with each support¬ 
ing sequence consisting of only one treatment. Note that ^ is the same 
Schur complement of as is that of V^. Dehne = L'^VsLi and 
Qs,i = LqVsLq. Let i/g,i(x) be the same function of Vgq as q^{x) is that of 
V^, and qs,i{x) be the same function of as qs{x) is that of Vs- Similar 
arguments for Theorem 1 yield the following theorem. 

Theorem 2. In estimating 9 under model (2), a pseudo symmetric mea¬ 
sure f with £'qV^£q > 0 is optimal in the following cases. In each case, the ps 
in ^ is positive only if s reaches the maximum therein. 

(i) If det{Q^^i) = 0, then ^ is Ag-optimal if and only if 

qi,i{x^,i)~'^qs,i{x^,i) + {t- 2){£'V^£)-‘^£'Vs£ _ ^ 
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//det(V^^i) > 0, then ^ is Ag-optimal if and only if 

r,,idet(Qg,i)/det(yg,i) + {t - 
s&s det(Q^^i)/det(V"^^i) + {t — 2){i'V^£)~^ 

(ii) //det(Qg^i) = 0, then ^ is Vg-optimal if and only if 

( 1 , t-2i'Vs£\_ 

^is\t-lq^^i{x^^i) t-li'V^iJ 
//det(V^^i) > 0, then ^ is Vg-optimal if and only if 

( Tsi t-2i'Vs£\ 

(iii) If det{Q^^i) = 0, then ^ is Sg-optimal if and only if 

max-—^—- = 1. 

6G5 

//det(V^^i) > 0, then ^ is Sg-optimal if and only if 

maxr<j i = 1. 
sg<S ’ 


= 1. 


(iv) If det{Q^^i) = 0, then ^ is Tg-optimal if and only if 

qs,i{x^,i) + {t-2)eVsi 
q^^i{x^^i)-\-{t - 2)i'V^i 

//det(V^^i) > 0, then is Sg-optimal if and only if 

Vs,I det(y^4)/det((5g4) + (t - 2)SVsl _ 
dei{V^,i)/dei{Q^,i) + {t-2)SV^l ~ ' 

Here and = -SV^Iq! 


Remark 5. With arguments similar to those in Remark 1, Theorem 2 
gives a comprehensive list of conditions to judge the optimality of pseudo 
symmetric measures with > 0 for estimating 6. Equivalence conditions 

for the four criteria could be easily derived when = 0, where we need 

to consider whether the cases of are equal to 0 or not, separately. We 

omit the details due to limit of space. 


Remark 6. If the within-block covariance matrix S is of type-H, = 

0 implies = 0, and thus = 0 in view of equation (22) below. Similarly, 
= 0 also results in = 0. Therefore, except for measures with each 
supporting sequence consisting of only one treatment, we have > 0 

and > 0 for any type-H matrix S. 
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The following proposition shows that the values of Ai and A 2 are irrele¬ 
vant to the determination of the Tg-optimal pseudo symmetric measures for 
estimating 6. 

Proposition 4. For any measure, we have 

qli = (1 + Ai + A2)^min[(l,x,y)rV^r(l,x,y)'], 

^iV 

where 



\0 0 1 / 


Proof. Let 

/I + A 1 + A 2 0 0 \ 

A = (1-|-Ai-l-A 2 ) ^ I Ai 1-|-A 2 —Ai j . 

\ A2 —A2 1 + Ai y 

Note that Li = (£,^ 05 ^ 1 ) = (1 + Ai -|- A 2 )rA. From Proposition 3 in Kunert 
and Martin (2000), we have 

=min[(l,x,y)P^,i(l,x,y)'] 

= (1 -L Ai -L A 2 )^ min[(l, x, y)A'rVgrA(l, x, y)'] 

= (1 -L Ai + A2)^min[(l,x,?/)r'Pgr(l,3:,?/)']. 
x,y 

The last equality uses the fact that for all possible values of x and y, 
A(l,x,y)' and (l,x,y)' share the same vector space. □ 

4. Optimal designs for the undirectional model. In many applications, 
it is reasonable to assume Ai = A 2 := A; that is, the neighbor effects do not 
depend on whether they are from the left or right. See Draper and Guttman 
(1980), Besag and Kempton (1986) and Filipiak (2012), for examples. Under 
this condition, model (2) reduces to 

(19) Irf = InkF + U/3 + {Td + ALd -|- XRd)T + e. 

The information matrix for r under model (19) is 
Cd{T) = 

The information matrix of a measure ^ for r is C^{t) = n~^Cd{T)- Also we 
consider only optimal measures in the pseudo symmetric format. 

Define I 2 = (0,1,1)' and L 2 = {^,^ 2 ), where i is defined in Section 2 with 
the value of (1, A, A)' here. Let U ^,2 = <3^,2 = Us ,2 = A'aU.As, 

Qs,2 = ^ 2 ^s ^2 and ?! 2 ~ Similar to Proposition 3, we 

have the following. 
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Proposition 5. For a pseudo symmetric measure the information 
matrix C^{t) has eigenvalues of Q, {t — and {t — with 

multiplicities of 1, 1 and t — 2, respectively. Moreover we have q '^2 — 

Note that if = 0, then q^ 2 ~ = c^oo hence C^{t) = c^ooBt/{t — 

1). By arguments similar to those in Theorem 1, we obtain the following re¬ 
sult. 


Theorem 3. In estimating r under model (19), a pseudo symmetric 
measure f, is optimal in the following cases. In each case, the ps in ^ is 
positive only if s reaches the maximum therein. 


(i) If ( 5^,2 = 0, then ^ is universally optimal if and only if 

CsOO 


max ■ 


= 1 . 


sG<S C^OO 

(ii) If det(V^^ 2 ) > 0, then is Ag-optimal if and only if 


max 

s^S 


r,,2QG2/detCPg,2) + {t - 2){i%I)-H'Vsi ^ 


Q^,2/detiV^,2) + {t-2){FV^i)-^ 


^ is Vg-optimal if and only if 


max 


r,,2 ^t-2£'Vsi 


s£S \t — 1 t — 1 i'V^i 


= 1 . 


^ is Sg-optimal if and only if 


maxTs 2 = 1- 
sG<S ’ 

^ is Tg-optimal if and only if 

rs,2det{V^^2)/Qi,2 + {t-2)eVs£ 

“s' det(y5,2)M,2 + (i - 

(hi) Otherwise, ^ is not optimal. 

Here ^ 5,2 = tr( 14 , 2 P^ 72 ^) - Qs, 2 Q^^l- 

It is easy to verify that for any measure, 

(20) = min[(l,a;)^^, 2 (l,a;)'] = min[(l,x,x)I7(l,x,x)']. 

Therefore, the value of A is irrelevant to the search of T^-optimal pseudo 
symmetric measures for estimating r. 

Next, we consider the total treatment effect 6. With the reason we ex¬ 
plained earlier, we shall assume 1 -|- 2A 7 ^ 0. The information matrix for 6 
under model (19) is 

CdiO) = (1 + 2A)-2M7;,^;,pr^(LI|M2,_i,_i0)Mi,yA. 
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For a measure its information matrix for 6 is given by C^{0) = n ^Cd{G)- 
Now we define £3 = (2, —1, —1)' and L 3 = (■^,■^ 3 ). Let = L'^V^L^, = 

ys ,3 = L'^VsLs, Q ,,3 = W 3 and Also 

we have the following. 


Proposition 6. For a pseudo symmetric measure the information 
matrix C^{9) has eigenvalues of 0, (l + 2A)“^(t —and (l + 2A)“^(t — 
with multiplicities ofl, 1 and t — 2, respectively. Moreover we have 

Qh < 


Note that if = 0, then 3 = i'V^i = (1 + 2A)^c^oo and hence C^{9) = 
c^ooBf/it — !)• Similar to Theorem 3, we obtain the following theorem. 


Theorem 4. In estimating 9 under model (19), a pseudo symmetric 
measure ^ is optimal in the following cases. In each case, the ps in ^ is 
positive only if s reaches the maximum therein. 


(i) If < 5^,3 = 0, then ^ is universally optimal if and only if 


ses C500 

(ii) //det(V^^ 3 ) > 0 , then f is Ag-optimal if and only if 

^ ,, rs,3/deijV^,3) + {t - 

™ Q^3/det(P^,3) + (t-2)(PP5^)-i 

^ is Vg-optimal if and only if 


max 


rs,3 

t-l 


+ 


t-2£’Vsi\ 

t -1 ev^i J 


f is £g-optimal if and only if 


= 1 . 


maxTs 3 = 1 . 

sdS ’ 

^ is Tg-optimal if and only if 

r,,3det(Pg,3)/Qg,3 + (t-2)fP,^ 
sg5 det(V^^3)/Q^^3 + (t — 2)PV^^ 

(hi) Otherwise, ^ is not optimal. 

Here rs ,3 = tr(P^, 3 P^“^) - Qs, 3 Q^^- 


It is easy to verify that for any measure, 

( 21 ) g |3 = min[(l,x)P^ 3 (l,a:)'] = (1 + 2 A)^ min[(l, x, x)rV^r(l, x, x)']. 

X X 
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Therefore, the value of A is also irrelevant in the search of £’g-optimal pseudo 
symmetric measures for estimating 0. 

Finally, we establish the connection between optimal measures for the 
directional and undirectional models if the within-block covariance matrix 
S is of type-H. 

Lemma 2 . If S is of type-H, we have q^2~ = 

(1 + Ai + A 2 ) 


Proof. Note that i? = if S is of type-H. For a sequence s = (ti,... ,ffc), 
dehne to = tk and tk+i = ti. Let kj be the frequency of treatment i appear¬ 
ing in s. Clearly, Y.\=ih = k. Let rus = k-^Y.\=i^h fs = Yli=i\=ti-i, 
9 s = Yli=i\=ti+i and hg = Yli=i\-i=ti+i ■ By straightforward calculations, 
we have c^oo — ^sii — Cs22 — k veig, c^oi — fs ’krig, C502 — 9 s 'krig and 
Csi2 = hg — mg. Since fg = gg, we have 

( k-nig fg-vUg 

fs i^s k Trig 

fs - mg hg- mg 

Note that I 4 = 0 if and only if s consists of only one treatment. From Propo¬ 
sition 3 in Kunert and Martin (2000), we have = minj;^y[(l, x, y)V^{l, x, y)']. 
Since (1,x,y)Vg(l,x, y)' is convex and exchangeable in x and y by equation 
(22) and = J2sesPs^s, h can achieve the minimum at some point of x = y. 
Therefore, = min 3 ;(l,x,x)Vg(l,x,x)' = q *^2 ™ view of equation (20). 

From Proposition 4, ^ = (1 -|-Ai -|-A 2 )^ mina;^y[(l,x,y)F'V^r(l,x,y)']. By 

equation (22) and Vg = J2s&sPs^s, we know (1, x, y)F'V^r(l, x, y)' is convex 
and exchangeable in x and y. Thus it can achieve the minimum at some 
point of X = y. Then (1 -|- Ai -|- A 2 )“^y|]^ = min,j,[(l!a;,a^)r^p 5 r(l!a^)^] = 
(1 -I- 2A)“^g|3 in view of equation (21). □ 

Theorem 5. If Ti is of type-H, a pseudo symmetrie measure is Eg- 
optimal for r (resp., 9) under model (19) if and only if it is £g-optimal 
for T (resp., 6) under model (2). Furthermore, if Ai = A 2 = A, the same 
result holds for Ag-, T>g- and Tg-optimal pseudo symmetric measures. 

This theorem is readily proved by using Lemma 2 and Propositions 2, 3, 

5 and 6 . 

5. Examples. For a sequence s = (ti,... ,tk), define the symmetric block 
of s as (s) = {as: a € O}. A symmetric block is an equivalence class, and 
hence S is partitioned into m -|- 1 symmetric blocks, say (sq), (si),..., (sm), 


fs - mg 
hg - mg 
k — mg 
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where Sj’s are the representative sequences in their own blocks. Without 
loss of generality, let (sq) be the symmetric block of sequences with identi¬ 
cal elements. For a measure ^ = {ps,s G 5 ), let = J2se{si)Ps = 

{P{si)^ ■ ■ ■ Since Vg is invariant for sequences in the same symmet¬ 

ric block, two pseudo symmetric measures with the same will share the 
same value of </>g(C)- By Remark 2 in Zheng (2013a), one can derive an exact 
optimal design in two steps: First, hnd the optimal P^, and then construct 
an exact pseudo symmetric design with that P^ by using some combinatory 
structures, such as type I orthogonal arrays [Rao (1961)]. See Aza’is, Bailey 
and Monod (1993) and Bailey and Druilhet (2014) for more techniques to 
construct exact pseudo symmetric designs. 

Note that Cgij =0, 0 <i,j <2, for any s G (sq). Given a measure ^ with 
Pl^so) > 0, one can always obtain a measure superior to ^ by replacing all 
sequences in (sq) with sequences not in the set. Therefore, the symmetric 
block (so) will be ignored in the following discusssion. 

In the sequel, we will determine the optimal P^ under model (2) through 
computer search based on Theorems 1 and 2. The one for the undirectional 
model (19) can be determined in a similar way by using Theorems 3, 4 and 5. 
The general algorithm for deriving the optimal P^ can be obtained by small 
modihcations of the algorithm in Zheng (2013b). For ease of illustration, we 
consider only 2 <t,k <5 and use the within-block covariance matrix to be of 
the form E = {Ii=j + /oli-j=±i(mod k))i<i,j<k- In the following examples, we 
take p in {0, —0.3,0.3}. Note that p = 0 implies E = that is, the errors are 
uncorrelated. First, let Ai and A 2 be nonnegative values from [0,1], and the 
negative case will be discussed later. All measures given below are pseudo 
symmetric measures. 

Cases k = 2 and 3. When k = 2, the symmetric block is (12). When k = 3, 
the symmetric block is (112) for t = 2, and those are (112) and (123) for t>3. 
By straightforward calculations, it can be verihed that the second smallest 
eigenvalues of C^{t) and C^(9) are both zero for any measure when k = 2 
and 3. Therefore, neither r nor 9 is estimable, and the optimal measures do 
not exist in these cases. This phenomenon is also observed by Bailey and 
Druilhet (2004) and Druilhet and Tinsson (2012) for the nonproportional 
interference model. 

Case of k = 4:. When t = 2, the four criteria become the same one. From 
Propositions 2, 3 and 4, it is known that the optimality of a measure for 
r or 0 does not depend on the values of Ai and A 2 . For r, we find that 
the measure with P(ii 22 ) = 1 is optimal for the three values of p. Next, 
consider the optimal measures for 0. If p = 0, the measure with P(ii 22 ) = 2/3 
and P(i 2 i 2 ) = 1/3 is optimal. An exact pseudo symmetric design with three 
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blocks based on it is given by 

/l 1 2 2 \ 

1 1 2 2 . 

\1 2 1 2 / 

In order to implement this design in practice, one has to adopt the ran¬ 
domization procedure as suggested by Azais, Bailey and Monod (1993). If 
p = —0.3, the measure with P{n 22 ) = 0-61 and P(i 2 i 2 ) = 0.39 is optimal. If 
p = 0.3, the measure with P(ii 22 ) = 0.76 and P(i 2 i 2 ) = 0.24 is optimal. 

When t = 3, the measure with P(ii 23 ) = 1 is optimal for r under the four 
criteria given all the values of Ai, A 2 and p. Consider the optimal measures 
for 6. If p = 0, the Ag-, Vg- and 7^-optimal measures vary for different 
values of Ai and A 2 . For all of them, there are two supporting symmetric 
blocks, that is, (1123) and (1213). Meanwhile, the former symmetric block 
dominates. The measure with ^( 1123 ) =2/3 and P(i 2 i 3 ) =1/3 is Tg-optimal. 
If p = —0.3 and 0.3, we observe that the supporting symmetric blocks are 
the same as those for p = 0 , except for the proportions of the supporting 
symmetric blocks. 

When t = 4 and 5, we find that the measure with ^( 1234 ) = 1 is optimal 
for both r and 9 under the four criteria, given all the values of Ai, A 2 and 
P- 

Case of k = 5. When t = 2, for both r and 0 we have the following. The 
measure with ^( 11122 ) = 0-8 and ^( 11212 ) = 0-2 is optimal for p = 0, the mea¬ 
sure with P(iii 22 ) = 0-71 and ^( 11212 ) = 0-29 is optimal for p = —0.3 and the 
measure with ^( 11122 ) = 0.90 and ^( 11212 ) = 0.10 is optimal for p = 0.3. 

When t = 3, first consider optimal measures for r. If p = 0, the optimal 
measures vary for different values of Ai and A 2 while the supporting sym¬ 
metric blocks are always (11223) and (12123). The proportion of (11223) 
is almost one for Ag-, Vg- and 7^-optimal measures and is 0.90 for the £g- 
optimal measure. If p = —0.3, the supporting symmetric blocks remain the 
same as those for p = 0 and (11223) still dominates. If p = 0.3, the measure 
with P(ii 223 ) = 1 is optimal under the four criteria for Ai, A 2 G [ 0 , 1 ]. For 6, 
we have observations similar to those for r. 

When t = 4, the supporting symmetric blocks are (11234) and (11223). 
When t = 5, the supporting symmetric blocks are (11234), (11223) and 
(12345). The optimal proportions and the dominating block may change 
for different values of Ai, A 2 and p. 

From Theorems 1 and 2, it is seen that as t increases, the equivalent 
conditions for optimal measures under Ag, Vg and Tg criteria tend to agree 
with each other. For example, take k = 5, p = 0, Ai = 0.1 and A 2 = 0.2. 
The measure with P(i 2345 ) = 1 is optimal under the three criteria for both 
r and 9 when t > 12. Meanwhile, the measure with P(ii 234 ) = 0.955 and 
P 1 J 12345 ) = 0.045 is Tg-optimal for both r and 9 as long as t > 5. 
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Table 1 

Efficiencies of optimal measures for r at 
(fc,t,p,Ai,A2) = (5,3,0,0.1,0.2) 


P<11223> 

A 

D 

E 

T 

0.98 

1 

0.99997 

0.98817 

0.99988 

0.99 

0.99998 

1 

0.98670 

0.99996 

0.90 

0.99265 

0.99213 

1 

0.99156 

1 

0.99988 

0.99997 

0.98496 

1 


Though the four criteria do not lead to the same optimal measure in 
general, the optimal measure under one criterion is typically highly efficient 
under the other three. Here the efficiency of a measure under a criterion 
is defined as the ratio of 4>g{^) to the maximum value among all measures. 
For the case of k = 5, t = 3, p = 0, Xi = 0.1 and A2 = 0.2, the efficiencies of 
optimal measures for r are shown in Table 1 and those for 9 under the four 
criteria are shown in Table 2. They all have efficiencies higher than 0.97. 
Furthermore, from the two tables we observe that the optimal measures for 
r also have high efficiencies in estimating 9 since they are almost the same 
as those for 9. 

From a practical viewpoint, the optimal proportions are sometimes too 
harsh for deriving exact designs. However, since the four criterion functions 
are continuous in the proportions, we could get a measure with good propor¬ 
tions in the neighborhood of the optimal one at the cost of a little efficiency. 
For example, when k = t = 5, p = 0, Ai = 0.1 and A2 = 0.2, the ^^-optimal 
measure for r is given by ^(11223) = 0.06 and ^(12345) = 0.94, which requires 
n to be a multiple of 50 at least. By rounding the proportions, we obtain a 
measure with ^(12345) = 1; which has efficiency higher than 0.99. An exact 
pseudo symmetric design with four blocks based on it is given by Azais, 


Table 2 

Efficiencies of optimal measures for 9 at 
(fc,t,p,Ai,A2) = (5,3,0,0.1,0.2) 


P<11223> 

A 

D 

E 

T 

0.93 

1 

0.99676 

0.98828 

0.98702 

0.99 

0.99556 

1 

0.98671 

0.99787 

0.90 

0.99859 

0.99213 

1 

0.97925 

1 

0.99307 

0.99981 

0.98215 

1 








20 


K. LI, W. ZHENG AND M. AI 


Bailey and Monod (1993) as 

/I 2 3 4 5\ 

1 3 5 2 4 

1 4 2 5 3- 

Vl 5 4 3 2/ 

In some occasions, the values of Ai and A2 could be negative. For ex¬ 
ample, a good fertilizer will possibly make a plant grow well so that the 
plant will compete with its neighbors for the sunlight, water and other re¬ 
sources in the soil. Suppose Ai,A2 € [—1,0). The Ag-, T>g- and 7^-optimal 
measures found by the computer program are different from the preceding 
ones for Ai,A2 G [0,1] in some cases. In estimating both r and 0, we ob¬ 
serve the following. When {k,t) = (4,3), the supporting symmetric blocks 
are (1123) and (1213) for Ag- and Pg-optimal measures, and are (1212) and 
(1213) for 7^-optimal measures. Contrarily, for Ai,A2 G [0,1], there is only 
one supporting symmetric block (1123) for optimal measures in estimating 
r, and the 7^-optimal measure in estimating 9 has two supporting symmetric 
blocks as (1123) and (1213). When {k,t) = (4,4), the ylg-optimal measure 
is still given by ^(1234) = 1- The supporting symmetric blocks are (1212), 
(1213) and (1234) for Pg-optimal measures, and are (1212) and (1213) for 
7^-optimal measures. But for Ai,A 2 G [0,1], there is only one supporting 
symmetric block (1234) for optimal measures under the three criteria. The 
details for other combinations of parameters are omitted due to the limit of 
space. 

6. Discussions. In this article, two proportional interference models are 
considered, in which the neighbor effects of a treatment are proportional to 
its direct effect. We investigate the optimal circular designs for the direct and 
total treatment effects. Kiefer’s equivalence theorems with respect to A, D, E 
and T criteria are established, based on which the search of optimal designs 
is easy to perform. Moreover, the connection between optimal designs for 
the two models is built. Examples are given to illustrate these theorems for 
several combinations of parameters. 

We now remark on directions for future work. Note that the number of 
distinct symmetric blocks will increase, at least geometrically, as the block 
size k grows. In such circumstance, it is unlikely that we could find the 
optimal proportions within a reasonable amount of time by using the cur¬ 
rent algorithm. Therefore, determining the forms of supporting symmetric 
blocks theoretically is vital to solving this problem. As a design theorist, the 
ultimate goal is to provide efficient or even optimal exact designs for any 
number of blocks. One way to achieve this is to further explore the construc¬ 
tions of exact pseudo symmetric designs. The other is to develop methods 
to build up efficient exact designs by modifying existing designs of smaller 
or larger size. 
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